o 
o 

(N 
bJO' 

; Quartz Tuning Fork: Thermometer, Pressure- and 

Viscometer for Helium Liquids* 



R. Blaauwgeers 1 , M. Blazkova 2 , M. Clovecko 3 , V.B. Eltsov 1 ' 41 ", 
O ' R. de Graaf 1 , J. Hosio 1 , M. Krusius 1 , D. Schmoranzer 5 , 

V ! W. Schoepe 6 , L. Skrbek 25 , P. Skyba 3 , R.E. Solntsev 1 , 

and D.E. Zmeev 4 

1 Low Temperature Laboratory, Helsinki University of Technology, Espoo, Finland 
2 Institute of Physics ASCR, Na Slovance 2, 182 21 Prague, Czech Republic 
3 Centre of Low Temperature Physics, Institute of Experimental Physics SAV, 

Kosice, Slovakia 

^Kapitza Institute for Physical Problems, Kosygina 2, Moscow 119334, Russia 
5 Faculty of Mathematics and Physics, Charles University, 
• Ke Karlovu 3, 121 16 Prague, Czech Republic 

6 Fakultdt fur Physik, Universitat Regensburg, D- 93040 Regensburg, Germany 

y—i '. 
> t 

Commercial quartz oscillators of the tuning-fork type with a resonant fre- 
£f~) | quency of ~ 32 kHz have been investigated in helium liquids. The oscillators 

OO ■ are found to have at best Q values in the range 10 5 — 10 6 , when measured 

o 

q | temperature operation the sensor has to be preselected. We explore their 

properties in the regime of linear viscous hydrodynamic response in normal 
and superfluid ^He and 4 He, by comparing measurements to the hydrody- 
namic model of the sensor. 

; PACS numbers: 47.80.-hv, 67.90.-hz, 85.50.Ly 

q ' KEY WORDS: quartz tuning fork, thermometer, viscometer, pressure sen- 
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*U 1 1. INTRODUCTION 

Quartz tuning forks are commercially produced piezoelectric oscillators 
meant to be used as frequency standards in watches. An extensive litera- 
ture describes their use for a large number of other additional applications.^ 



* Dedicated to former JLTP editor Frank Pobell 
^E-mail: ve@boojum.hut.fi 
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They have also been employed in liquid He temperature measurements P This 
study was inspired by the expectation that industrially produced quartz os- 
cillators, with a calibrated standard frequency of 2 15 Hz (= 32 768 Hz) at 
room temperature, would be reasonably identical and could be used as sec- 
ondary thermometers without need for recalibration. 

We have performed measurements on four different forks of identical 
dimensions, but produced by different manufacturers. It turns out that 
without preselection at LHe temperatures the results cannot be reduced 
on a common temperature dependence. For instance, the resonance width 
A/ vac measured in vacuum below 1.5 K proved to be 0.06 Hz, 0.5 Hz, 1.4 Hz, 
and 0.06 Hz for these four sensors. This measure of the intrinsic dissipation, 
which limits the response of the device to its environment at the lowest 
temperatures, cannot be determined from room temperature measurements. 
It remains to be seen if simple preselection criteria can be worked out, to 
narrow down the variation in oscillator characteristics. 

However, the quartz tuning fork offers other important advantag es as a 
sensor of its cryogenic environment. Forks are cheap and readily available, 
they are robust and as such easy to install and to use, they operate at a 
higher frequency than most other vibrating sensors, and they are highly sen- 
sitive indicators of the physical properties of the medium in which they are 
immersed. Thus they provide handy in situ information about the condi- 
tions in a sample container at the far end of a low-conductance filling line, 
which is helpful during flushing, filling, emptying, and in general, for repro- 
ducible monitoring of pressure and temperature changes. A major advantage 
in many applications is that to drive these piezoelectric devices no magnetic 
fields are needed and that they, in fact, are highly insensitive to themPSl 

In this report we explore the linear response of the quartz tuning fork 
in He liquids at low excitation, in the regime of viscous hydrodynamics. The 
purpose is to compare the measurements to a hydrodynamic model which 
could explain the measured results. For this Sec. 2 studies the oscillator 
properties of the tuning fork in vacuum. In Sec. 3 the influence of the 
surrounding medium is incorporated, Sec. 4 discusses briefly the practical 
measurement, and later sections describe the measured results in vacuum as 
well as in 3 He and 4 He liquids, by comparing the data to the physical model. 
We postpone to a later occasion the analysis of nonlinear effects and of 
the lowest temperatures with collisionless motion of excitations. This latter 
aspect, the creation and detection of excitations and of quantized vortices 
in the T — > temperature limit, is of great current interest. A large amount 
of new information has been discovered on vortex properties using vibrating 
resonators: (i) spherespl grids,^ and wires^ in 4 He-II, (ii) grids 7 and wires^ 
in 3 He-B, and (iii) wires^ in 3 He- 4 He mixtures. One might hope that the 
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quartz tuning fork could be used for similar measurements. 



2. CHARACTERISTICS OF THE VIBRATING TUNING 

FORK 



At sufficiently small oscillation amplitudes the fork can be described as 
a harmonic oscillator subject to a harmonic driving force F = Fq cos(u>i) and 
a drag force with linear dependence on velocity. The equation of motion is 
given by 



We have here four parameters, namely the effective mass m (of one leg), the 
drag coefficient 7, the spring constant k, and the amplitude of the driving 
force Fq. The effective mass and the drag coefficient depend on the medium 
around the oscillating fork. The solution of this differential equation is well 
known: It can be written as x(t) = x a (uj) sin(ut) + Xd(u>) cos(wt), where x a 
and Xd are the absorption and dispersion, respectively. The mean absorbed 
power (F dx/dt) = Fqux^/2 is at maximum at the resonant frequency 



It is convenient to introduce the quality factor 




as the ratio of the resonant frequency ujq to the frequency width Auj = 7, 
where Auj is the full width of the resonance curve at half of the maximum 
power. 

The geometry of the fork is sketched in Fig. ^ It is characterized by the 
length C, width W, and thickness T of a leg. The relevant vibration mode 
is the basic antisymmetric mode, i.e. the one where the two legs of the fork 
move in antiphase along the direction of T. Taking the known elasticity 
modulus E of quartz, E = 7.87 • 10 10 N/m 2 , the spring constant is given 



2.1. Mechanical Properties 
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and the effective mass of one leg in vacuum is 



m vac = 0.24267 pqvCWT, 



(2) 
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Fig. 1. Sketch of the quartz tuning fork. 



where we use for the density p q = 2659 kg/m 3 , the density of quartz (ne- 
glecting the electrodes on the legs). 

We use forks with C = 3.12 mm, W = 0.352 mm, T = 0.402 mm, and 
T> = 1.0 mm. For this geometry we find k = 1.48 • 10 4 N/m and m vac = 
2.85 • 10~ 7 kg. This gives / = ujq/2h = 36293 Hz, which is 11% larger than 
the manufactured value of 32768 Hz at room temperature. The discrepancy 
with the theoretical expression is most likely due to additional weight of 
the evaporated electrodes, dependence of the elasticity modulus of quartz 
on the orientation with respect to the crystallographic axes, and deviations 
in geometry between the real fork and the model. At room temperature 
in vacuum the measured width of the absorption curve at half power is 
typically A/ ~ 1.2 Hz, which gives Q = /q/A/ ~ 2.7 ■ 10 4 . In room air the 
width increases to A/ ~ 4.3 Hz. With decreasing temperature the resonant 
frequency diminishes and the Q value increases. Our vacuum measurements 
at LHe temperatures will be presented in Sec. ED 



2.2. Electrical Properties and Calibration 

The fork is excited with ac voltage U = Uq cos(u;i) while the frequency 
is slowly swept through resonance. The signal received from the fork is a 
current I owing to the piezoelectric effect. The stresses due to fork deflection 
induce charges and thus the current is proportional to the derivative of the 
fork deflection, i.e. to the velocity: 



(3) 



Quartz Tuning Fork 



where a is the fork constant. Its theoretical value is given by^ 

a = 3d 11 E(TW/£) , (4) 

where d\\ = 2.31 • 10 _12 m/V is the longitudinal piezoelectric modulus of 
quartz. For our forks a = 2.47 • 10~ 5 C/m. The oscillation amplitude is 
therefore known theoretically, but in practice the fork is usually calibrated 
optically with interferometric techniques! 10 * 11 ^ Typically, only ca. 30% of 
the theoretical current sensitivity is achieved. In a cryogenic setup, where 
the fork is mounted inside a sample container in the heart of the cryostat, 
a direct measurement of the fork constant a is complicated. Instead, we 
prefer to calibrate our forks by comparing the mechanical oscillator with 
the equivalent electrical RLC series resonance circuitP The corresponding 
differential equation for the current is 

d 2 I Rdl_ _Z_ ldU 

Comparing Eqs. (J5J) and Q we see that ujI = l/(LC), 7 = R/L, and 
using Eq. (J3J), 1/L = (Fq/Uq) a/m. Additionally we have the condition 
that the dissipated power at resonance has to be equal for both equations: 
The electrical power Uq/(2R) drives two legs of the fork which dissipate 
2 • F 2 /(2m7). Thus we have a closed set of equations which allows us to 
connect the electrical and mechanical properties of the fork via the fork 
constant a: 

F = (a/2) U , (6) 
R = 2m 7 /a 2 , (7) 
L = 2m/a 2 , (8) 
C = a 2 /(2k) . (9) 

Experimentally the fork constant a can be determined using Eq. Ij7j). which 
can be rewritten as 

^ . do) 

Here Alo is determined from the width of the resonance curve while 1/R is 
the linear slope of the experimental Io{Uq) dependence, where Iq is the cur- 
rent amplitude at resonance. The only parameter which cannot be directly 
determined from the experiment is the effective mass m. However the theo- 
retical value of the effective mass (Eq. EJ) seems to be fairly reliable because 
of the close agreement between theoretical and experimental values of the 
resonant frequency. The example of our fork response in vacuum in Sec. ED 
leads to a = 8.13 ■ 10~ 6 C/m, which amounts to 33 % of the theoretical value. 
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In our experiment the absorption I a (k>) and dispersion Id{uo) compo- 
nents of the current I(t) = I a cos(wi) + Id sin(wt) are measured separately 
with a lock-in amplifier (see Sec. El for details). The theoretical resonance 
curves, 

a 2 Uo mjuj 2 
a = ~T~ (m 7 cu) 2 + H 2 - k) 2 ~ {Auo) 2 uj 2 + (uj 2 - a;2)2 ' 

a 2 Uo uj{muj 2 — k) IqAuj uj (uj 2 — ujq) 

d = ~2~ (mjuj) 2 + (muj 2 - k) 2 = (Aw) 2 w 2 + {uj 2 - uj 2 ) 2 ' ( ' 

can be fit to the experimental response to determine the parameters which 
enter Eq. (|T|). In particular, the absorption component I & {to) reaches its 
maximum value Iq at a frequency which is exactly ujq and the full width Auj 
of the absorption curve at 1 /2 of the maximum height Jo gives exactly 7. If 
the fork constant a is known, then m = a 2 Uo/(2IoA.uj) and k = mojQ can be 
determined independently. 



I (Acu) 2 lu 2 



3. INFLUENCE OF SURROUNDING MEDIUM ON THE 

OSCILLATING FORK 

3.1. Hydrodynamic properties 

In this section we outline the basic properties of the oscillatory boundary 
layer flows, to understand how the vibrating fork works as a detector. The 
classical viscous flow around a submerged oscillating body^l is rotational 
within a certain layer adjacent to the body, while at larger distances it 
rapidly changes to potential flow (if there is no free liquid surface or solid 
surface in the vicinity of the oscillating body). The depth of penetration of 
the rotational flow is of order 




where uj is the angular frequency of oscillation while 77 and v = rj/p are the 
dynamic and kinematic viscosities of the fluid with density p. 

As a result of the oscillatory motion of the body through the liquid, the 
body experiences a force which has components proportional to the velocity 
of the body v (drag) and to its acceleration v (mass enhancement): 

T = bv + fhv. (14) 

To determine the values of b and fh generally a full solution of the flow 
field around the oscillating body is required. Simplifications are possible in 
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two limiting cases, which depend on the relative magnitudes of the charac- 
teristic size of the oscillating body £, oscillation amplitude xo, and viscous 
penetration depth 5: 

1) £ -C 5 and loxq£/v <C 1: In this case the flow at any given instant 
can be regarded as steady - as if the body were moving uniformly with its 
instantaneous velocity. As a rule, this case does not apply to the oscillating 
cryogenic flows considered here. 

2) £ S> 5 and £ 3> xq: In this case the layer of rotational flow around 
the body is very thin while in the rest of the fluid the flow is potential. This 
case is directly applicable to quartz tuning forks in 3 He and 4 He liquids: 
The kinematic viscosity of normal liquid 4 He-I above the A-point is ~ 
2-10 -4 cm 2 /sand of normal 3 He above the superfluid transition v% ~ lcm 2 /s. 
At 32 kHz we get the penetration depths £4 as 0.4 pm and £3 ~ 30 /im while 
£ ~ 400 /im for our forks. Moreover, the oscillation amplitude would reach 
the leg thickness T only at a very high velocity of order lm/s (c/. Fig. EJ). 
Note that other oscillating objects such as spheres or wires may not always 
be in this flow regime since they usually have smaller characteristic size and 
smaller oscillation frequency. 

When the conditions for case (2) above are valid, a major contribution 
to both b and m in Eq. (|14[) is found by solving the potential flow field u 
around the bodyP^ In particular, b is expressed as 

-ffics, (15) 

where vq and uq are the amplitudes of the velocities of the body and the 
flow, the integral is taken over the surface of the oscillating body, S is the 
surface area of the body, and C is some numerical constant which depends 
on the exact geometry of the body. For example for a sphere C = 3/2, while 
for an infinitely long cylinder oscillating perpendicular to its axis C = 2. 

The largest contribution to mass enhancement rh comes from the po- 
tential flow around the body and can be expressed through the mass pV 
of the liquid displaced by the body of volume V. A smaller contribution 
is caused by the fact that the viscous drag force experienced by the body 
is usually phase shifted with respect to the velocity of the body. This can 
be interpreted such that a volume of order SS of the liquid is clamped to 
comotion with the oscillating body. Thus for fh we can write 

m = f3pV + BpS6 , (16) 

where (5 and B are again geometry-dependent coefficients. For example, for 
a sphere^ (5 = 1/2 and B = 3/4; for an infinitely long cylinder with elliptic 
cross section^ (5 = r±/r», where r± and r» are the lengths of the axes which 
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are perpendicular and parallel to the oscillation direction, respectively; for 
a rectangular bearrPl (3 = (7r/4)a±/a^ , where aj_ and ay are the lengths of 
the sides which are perpendicular and parallel to the oscillation direction, 
respectively. 

We are not aware of rigorous calculations of the parameters 0, B, and 
C for a tuning fork. A single oscillating beam has been considered before 
in great detail. — However, the presence of two legs in close vicinity of each 
other significantly affects the potential flow field and changes, for example, 
the (5 parameter P Thus we consider /?, B, and C as fitting parameters, to 
be determined for a particular fork from the experiment. This approach was 
previously used, for example, in Ref. ED f° r a vibrating reed in liquid 4 He 
and it seems to be provide the first step for understanding the experimental 
results. 



3.2. Hydrodynamic model of sensor 

The addition of the force T to the equation of motion of the fork Q 
leads to a reduction in the resonant frequency and an increase in the width 
of the resonance curve: 

w = w 0vac( m vac/m), (17) 
7 = 7vac(™<vac/"1<) + b/m, (18) 

where m = m vac + m is the effective mass of the oscillating body immersed 
in the fluid. For convenience we redefine the fork parameters (3 and B from 
Eq. (|16|) as relative to the effective fork mass in vacuum: 

m = m vac £-\p + B5{S/V)} . (19) 

Ignoring the vacuum resonance width A/ vac , we finally obtain the depen- 
dence of the resonant frequency fo and of the full width of the absorption 
curve at half height A/ on the fluid density and viscosity: 



/0vac ^ 1 + ^ (P + B ) , (20, 



fo J Pq V V V 7rp/o 



\2 



Af = - PnA C S Uo/fM . (21) 
2 V vr m vac 

Here V = TWC and S = 2(T + W)C. 

Eqs. Q2UJI and (|21|) can now be used to determine experimentally the 
hydrodynamic parameters from measurements in a fluid with known p and 
rj. Once the parameters are known, the fork can be used for measurements 
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of p and rj, in principle, for any other medium and thus as a pressure and 
temperature sensor, if p(P, T) and rj(P, T) are known. From this point of 
view, the width Af is especially useful for measurements as it requires cali- 
bration of only one parameter (C/m vac ). In Eq. (|2U|) the factor multiplying 
the parameter B is small and so, even here, often only one parameter {(5) is 
of major importance. Unfortunately, our measurements indicate that the pa- 
rameters vary from one fork to the next and calibrations need to be checked 
(see Sees. 6 and 7). The calibration should be re-examined even when per- 
forming measurements with the same fork in widely differing conditions (see 
Sec. 6). 

3.3. Beyond the model of viscous hydrodynamics 

With decreasing temperature the mean free path of excitations increases 
in both 4 He and 3 He superfluids and the hydrodynamic description ceases 
to be valid as the normal fluid penetration depth grows beyond all relevant 
length scales. In 4 He-II the crossover to the ballistic regime takes place be- 
low IK and in 3 He-B below 0.3 T c . At low temperatures the drag is caused 
by the scattering of the excitations from the oscillating body. As the excita- 
tion density decreases with decreasing temperature the drag coefficient also 
rapidly decreases: b oc T 4 for phonons in 4 He-II, while b cx exp(— A/k^T) 
for rotons and for quasiparticles in 3 He-B, where A is the relevant energy 
gap. The crossover from the hydrodynamic to the ballistic regime has been 
described for a vibrating wire in Refs. |16J17j and for a vibrating sphere in 
|18j . In this work we are not discussing the ballistic regime further. 

In everyday life a tuning fork is used to create sound in air. In the case 
of a quartz tuning fork in LHe we might wonder whether the compressibility 
and the losses from sound emission need to be taken into account. A quartz 
fork operates at higher frequency than vibrating wires, grids, spheres and 
most other oscillating bodies. Therefore sound emission might be sizable 
for the fork while it is negligible in these other cases. The power loss from 
acoustic emission reduces the Q value and contributes to the width of the 
resonance curve Au; ac = i? ac /m, where the average power loss ^R^Vq has 
been expressed in terms of the so-called radiation resistance i? ac . 

A realistic calculation of the acoustic emission from a tuning fork is 
complicated. Clubb et al)^ suggest a model of two infinite cylinders oscil- 
lating at 180° out-of-phase and give for their quadrupolar acoustic field the 
radiation resistance 

11616 c 4 v ; 

This expression gives a very small contribution to the resonance width and, 
owing to the high powers in which the different quantities appear, quantita- 
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Fig. 2. Circuit diagram for tuning fork measurements. 



tive comparison with experiment is so far inconclusive. In our measurements 
on normal 3 He a small temperature independent constant contribution to the 
resonance width is distinguishable at high temperatures (see Sec. 6). How- 
ever, its magnitude is not in agreement with Eq. (I22H and at this point the 
presence of the acoustic loss term remains unclear. 



4. SENSOR PREPARATION AND MEASUREMENT 

A commercial quartz tuning fork comes in a vacuum-tight sealed metal 
can which has to be partly or entirely removed, to probe the flow properties 
of the surrounding medium. The pair of leads for exciting and sensing the 
oscillator is magnetic. In magnetically sensitive applications they have to 
be removed and changed to nonmagnetic leads. For measuring the current 
in Eq. © one usually uses the current input of a phase-sensitive lock-in 
amplifier. The excitation voltage is supplied by a high resolution digital 
generator, which also provides the reference signal for the lock-in amplifier 
(see circuit diagram in Fig. |2J). 

In Fig. El the results are plotted from room temperature test measure- 
ments, to compare to Eqs. (|20(l and (|21|) . The measured nearly linear de- 
pendence of the resonant frequency and the square root dependence of the 
resonance width versus applied pressure P follow directly from these equa- 
tions, assuming that p oc P and r] does not appreciably change with P. In 
these measurements we also tested if the results showed variations depending 
on whether (i) only a small hole is made in the encapsulating can, or (ii) the 
entire top surface of the can is ground away, or (iii) the can is completely 
removed. No obvious qualitative differences were observed, which is as ex- 
pected, since for all data in Fig. |31 the penetration depth S is much smaller 
than the fork dimensions. 



Quartz Tuning Fork 




20 



32700 ■ °' 



He 



^ 32650 ■ 
a, 




- 15 



< 32600 ■ 




32550 : 







5 



10 



15 











20 



Pressure P, bar 



Fig. 3. Room temperature tests of sensor sensitivity. (Main panel) Reso- 
nance width A/, with the vacuum width A/ vac subtracted (right vertical 
scale), plotted versus applied pressure in gaseous nitrogen (T = 22.5°C) and 
helium (T = 23°C). The solid line is the fitted square root dependence, as 
expected from Eq. (|21|). (Inset) The corresponding resonant frequencies /o 
are almost linear with pressure, as expected from Eq. ()20[) , The data for N2 
and He have been measured with two different forks. 

In Fig. |3] a similar measurement has been performed in liquid 4 He at 
4.2 K as a function of pressure. Owing to the non-linear dependence of the 
density of liquid He on applied pressure the resonant frequency shows linear 
dependence only after converting pressures to densities. For these results 
thermal aging is of some concern. We tested the stability of the resonant 
frequency and width of five different forks by cycling them between 300 
and 77 K, with their cans completely removed. For a virgin sensor resonant 
frequency shifts of < 0.2 Hz and changes in width < 0.3 Hz are typical. For 
most applications such shifts are negligible, but if more stable reproducible 
results are required, then the sensor should be thermally cycled. After a few 
cycles the changes are considerably reduced, but only after 30 - 50 cycles the 
results become stable. It should be noted that larger changes can result from 
bending the leads of the fork or from making new solderings to the leads, 
presumably because new strains are imposed via the electrodes on the quartz 
surface. Such changes, which can be < 0.7 Hz in both resonant frequency 
and width, are generally larger than those caused by the removal of the can 
by grinding. This suggests that a fork, which is used as thermometer, should 
be handled with care. 
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Fig. 4. Resonant frequency /q in liquid 4 He at 4.2 K plotted versus applied 
pressure (bottom) and density (top). The tuning fork was in its original can, 
but with the fiat top surface of the can ground away and A/ vac = 0.06 Hz. 
The pressure was converted to density using the HEPAK packageP The 
solid line is a linear fit through the data. 



A further practical consideration in connection with Fig.|l]is that these 
results could only be obtained after the tuning fork was installed in an iso- 
lated sample container in controlled conditions. In a LHe bath in an open 
dewar surface conditions on the fork may change because of adsorbed gas 
or particles floating in the bath after a LHe transfer. Typically in such con- 
ditions the resonance width does not stay constant, but gradually increases 
during a long run. Occasionally small step changes in the resonant frequency 
are observed which could arise if air flakes stick on the fork, for instance. 
Thus for accurate and reproducible readings the fork should preferably not 
be used in technical helium. 

Fig. El shows two examples of the quartz tuning fork as a practical mon- 
itoring device of temperatures in a nuclear cooling cryostat! 21 ^ 21 * The two 
traces illustrate the range which the resonance width A/ traverses when the 
cryostat is taken through its cooling cycle. The bottom trace represents a 
sensor in the mixing chamber, while the top trace monitors one on the nu- 
clear cooling stage. The most prominent features in these traces are abrupt 
anomalies: In the bottom trace from disconnecting the superconducting heat 



Quartz Tuning Fork 




Time, min 



Fig. 5. Two quartz tuning forks monitoring temperatures in a nuclear cooling 
cryostat. The bottom trace shows the sensor inside the mixing chamber of 
the dilution refrigerator in the concentrated 3 He- 4 He solution. The top trace 
gives the simultaneously measured width of the sensor on the nuclear cooling 
stage inside the liquid 3 He sample container at zero liquid pressure. 



switch between the mixing chamber and the nuclear cooling stage. Here the 
mixing chamber first starts cooling, but then warms up to a new higher level 
when the 3 He flow rate in the dilution refrigerator circulation is reduced. 
Simultaneously the demagnetization of the nuclear cooling stage is started. 
This means that in the upper trace the temperature starts decreasing and 
A/ increasing. When the cool down passes through the superfluid transition 
temperature T c , the width suddenly dives in a rapid decrease (see Sec. 6). 
These recordings show that the quartz tuning fork provides rapid and sen- 
sitive confirmation of the actions which are performed on the cryostat. 



5. PROPERTIES OF TUNING FORKS IN VACUUM 

In its original package the quartz tuning fork comes inside a vacuum- 
tight can. When the fork is cooled in this can, the resonant frequency de- 
creases and the Q value increases. At LHe temperatures the reduction in 
resonant frequency from the room-temperature value is about 70 Hz and the 
Q value approaches 10 6 . If the can is removed the bare fork behaves in vac- 
uum in similar manner: The reduction in the resonant frequency remains 
the same while the Q value is typically lower than for the fork inside its 
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Fig. 6. Response of the quartz tuning fork in vacuum at a temperature of 
about IK. (Left) Resonance curves for absorption L A (o) and dispersion 
(□) components of the current measured at fixed excitation level Uq. The 
solid curves are fits to Eqs. Qll|) and Background contributions from a 

shunting capacitance and leads are subtracted. (Right) Dependence of the 
amplitude Iq of the current at resonance on the excitation voltage amplitude 
Uq. Conversions to velocity and driving force via Eqs. (j2J) and (JSJ) are also 
shown. 



original can. However, the Q value varies greatly from one fork to the next. 
Typical Q values range from 2 ■ 10 4 to 5 ■ 10 5 . A similar large scatter was 
observed in previous reports."™ The reasons for such variability between dif- 
ferent forks are not known. In principle differences may result from slight 
damage (when the can is removed, for instance) or from dirt accumulated 
on the legs. In practice of course, the resonance width is much larger in 
LHe than the vacuum width and thus these problems do not affect our later 
analysis. However, if one is interested in the properties at the lowest temper- 
atures in the ballistic regime, then special care has to be invested in selecting 
forks with the narrowest possible vacuum width. 

The vacuum response of a tuning fork at ~ 1 K is shown in Fig. |3 This 
fork was later used for measurements in 3 He liquid at zero pressure. The 
data in Fig. were measured at different excitation levels in the 3 He sample 
container while the 3 He pressure was less than 1 mbar at a temperature of 
around IK. As mentioned above (Sec. 2.2), from the resonance character- 
istics measured at various excitation voltage amplitudes Uq it is possible to 
determine the fork constant a using Eq. (|10|). In this case a = 8.13-10 -6 C/m. 
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Now the driving force -Fo and the velocity amplitude can be calculated from 
Eqs. © and ©, respectively. Thus the measurement of Iq versus XJq in the 
right panel of Fig. H3 can be converted to a dependence of the tuning fork 
velocity amplitude on the driving force. 

As seen in Fig. H2 in vacuum the fork responds linearly up to a driv- 
ing force of order 100 nN, which corresponds to displacements of a few jiva.. 
Above this limit the amplitude in the motion of the legs starts to be sufficient 
for the response to become nonlinear. The cause for the nonlinear behavior 
is deformation in the sensor material during large-amplitude oscillation. The 
oscillation amplitudes of ions around the minimum of their potential become 
large enough such that anharmonic terms in the potential energy introduce 
nonlinear restoring forces in the ion motion. The overall effect is the ap- 
pearance of a nonlinear restoring force in the equation of motion, Eq. 
In addition at low temperatures other sources contribute to nonlinearities in 
large amplitude oscillation, such as slow strain release from defects in the 
oscillator material. These have been extensively investigated, for instance 
with vibrating wire resonators in vacuumP3 The nonlinear drive regime of 
the quartz fork is not discussed in this report. 



6. TUNING FORK IN 3 He 

Our measurements in liquid 3 He have been performed in two different 
cryostats and with different sensorsP^ In the measurements at 29 bar pres- 
sure, the fork parameters were /o V ac = 32705.05 Hz and A/ vac = 0.06 Hz 
in vacuum at LHe temperatures. In this setup temperatures above T c are 
measured by a melting curve thermometer which is mounted on the nuclear 
cooling stage. Below T c , the temperature readings are determined from nu- 
clear magnetic resonance frequency shifts of the 3 He samplePS The NMR 
reading is preferred below T c because it is measured directly from the liquid 
in which the tuning fork is also immersed. This minimizes thermal gradients 
between the thermometer and the fork. 

In the zero pressure measurements the fork had /o V ac = 32707.4 Hz and 
A/vac = 1.41 Hz. In this setup the temperature is determined with pulsed 
NMR on Pt powder immersed in the liquid 3 He sample. The NMR signal 
amplitude is calibrated using the known value of the superfluid 3 He transition 
temperature T c . The superfluid transition is indicated by the fork reading 
or by two additional vibrating wire resonators in the 3 He cell which also are 
used for thermometry. 

The temperature dependence of the resonance width A/ at 29 bar pres- 
sure is presented in Fig. [7| The measurements were performed in the linear 
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Fig. 7. Resonance width of the quartz tuning fork in liquid 3 He at 29 bar 
pressure. The measured data are marked with circles. A sharp reduction in 
the width is observed on cooling below T c and an abrupt discontinuity at 
the AB transition. The solid line is the predicted behavioi^ of a vibrating 
wire resonator with the same density p q and the same vacuum resonant 
frequency /o V ac as the fork. To produce a good fit to the quartz fork data 
the wire diameter had to be fixed to 0.25 mm which is comparable to the 
dimensions of the legs of the fork. 

drive regime of the fork with the maximum current not exceeding 4 nA (and 
thus velocities not exceeding 0.5mm/s). The plot demonstrates that the 
width changes rapidly in the range 0.8-40mK and, once calibrated, can 
be used as a thermometer. The measurements were extended to much lower 
temperatures than shown in Fig. [JJ the lowest observed resonance width was 
Af = 0.5 Hz. However, at these temperatures we have no other thermometer 
in the 3 He sample container to calibrate the fork. 

Fig. UJ can be divided in three temperature regimes: (i) normal 3 He 
above T c , where the width rapidly increases with decreasing temperature 
due to the Fermi- liquid behavior of the viscosity, rj cx 1/T 2 , (ii) superfluid 
3 He-A and (iii) 3 He-B phases, where the width decreases with decreasing 
temperature mainly because of the decreasing normal-fluid density. Fig. [7J 
provides an interesting comparison of the fork oscillator with the vibrating 
wire resonator. In Ref. the available theoretical and experimental infor- 
mation on vibrating wires has been combined in a computer program which 
calculates the response of a vibrating wire loop in normal 3 He and in su- 
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Fig. 8. Closeup of resonance width in Fig. in the temperature region 
of 3 He-A. During cooling the A— >B transition is supercooled, while during 
warming the B— >A transition occurs close to the thermodynamic equilibrium 
temperature Tab of this first order phase transition. The discontinuous jump 
in resonance width at the transition is mainly caused by the change in p n . 
The solid lines depict the vibrating wire model of Fig. [7| 

perfluid 3 He-B for a resonator with known wire diameter, wire density, and 
resonant frequency in vacuum. In Fig. [7|the response of a fictitious wire loop 
with the density of quartz and the frequency of the tuning fork in vacuum 
has been fitted to the experimental data with the wire diameter as a fitting 
parameter. The fit is remarkably good and gives a reasonable value for the 
wire diameter which is of order of the thickness of the fork leg. The same 
fit does not reproduce exactly the temperature dependence of the resonant 
frequency owing to the difference in the (3 factor between the wire and the 
fork (c/. Eqs. (|20|) and (|21jl). The comparison is used here to emphasize 
that the fork thermometer is very comparable to the vibrating wire in this 
temperature range of liquid 3 He. 

Vibrating wires are not generally used as thermometers in 3 He-A owing 
to their texture-dependent non-reproducible response. Our limited expe- 
rience with forks shows that for a given fork the width in the A phase is 
reproducible, independently whether one enters the A phase from the nor- 
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Fig. 9. Resonance width at two different pressures in normal 3 He as a func- 
tion of inverse temperature. The experimental data is shown as circles. The 
solid lines are linear fits of Af versus 1/T. For comparison, the dashed 
curves show the effect from a viscosity anomaly close to T c which was mea- 
sured with a vibrating wire resonator in Ref. |25| . 



mal phase or from the B phase direction, as seen in Fig. |HJ Probably the 
larger dimensions of the fork legs fix the orientations in the order parameter 
texture such that the response becomes reproducible if the oscillation am- 
plitude remains small compared to the leg dimensions. Nevertheless, more 
measurements are required to establish whether tuning forks can be used as 
accurate secondary thermometers also in 3 He-A. 

In Fig. Elthe fork properties are analyzed in more detail in normal 3 He. 
Since the viscosity of normal 3 He varies as r\ oc T~ 2 , we expect the width 
to depend on temperature as Af oc T -1 , Eq. (|2"Tj) . This dependence is 
indeed observed in the experiment. However, a constant addition to the 
resonance width is also present in the experimental data (seen as a non-zero 
intercept on the vertical scale in the two panels of Fig. EJ). This additional 
temperature-independent contribution to the width is 21.2 Hz at 29 bar and 
3.2 Hz at zero pressure. In Fig. this effect appears as a tendency towards 
a constant width at high temperatures. 

The origin of the temperature-independent contribution to the width 
in normal 3 He is not clear. Clubb et alP observed the same effect in 3 He- 
4 He mixtures and attributed it to acoustic emission. Their model, Eq. I|22j>. 
gives an orders of magnitude smaller value than the measured one in our 
case. Moreover, irrespectively of the model one would expect that the losses 
from sound emission will be smaller at high pressures (since the sound ve- 
locity increases with pressure), while the two measurements in Fig. ^Jshow 
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Fig. 10. Resonance frequencies of the two quartz tuning forks from Fig. |5] 
in normal 3 He, plotted as [(/ovac//o) 2 — IJPq/p versus l/T. The symbols 
represent experimental data. The solid lines are fits to Eq. (|'20j) . assuming 
7] cx T~ 2 . 

opposite behavior. Since these have been performed with different sensors, 
the possibility remains at this point that the main part of the temperature 
independent width in normal liquid 3 He depends on the fork with which it 
is measured. 

Several reports on vibrating wire measurements mention a viscosity 
anomaly in normal 3 He: an unexpected reduction in the viscosity close to 
T c from the T~ 2 behaviorP3 Our data show no sign of this anomaly: The 
width Af changes exactly proportional to T _1 until T c (Fig.EiJ). The reason 
for this difference is not clear. Possibly the fork owing to larger size and 
higher frequency operates in a different hydrodynamic regime than typical 
vibrating wire resonators. In particular close to T c , where the viscosity of 
3 He is the highest, the viscous penetration depth becomes comparable to 
the characteristic size of the oscillating object, especially at low pressure. 
For example for a fork at P = and T = 1 mK, 5 is about the inter-leg 
distance T>. Thus an interpretation of the results in terms of the simple 
model presented in Sec. 3 may not be justified. The question which kind of 
viscometer is more appropriate for 3 He at temperatures close to T c requires 
further analysis. 

From the linear fit in Fig. |2]the fork parameter C in Eq. (|21[) can be de- 
termined. Using viscosity data from Ref. [221 (omitting the viscosity anomaly 
close to T c ) we get C = 0.57 for the fork used at zero pressure and C = 0.64 
for the fork used at 29 bar . 

The relative change in the density of liquid 3 He between and 40 mK 
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is about 10 -4 , according to Ref. |28| . Thus the largest contribution to the 
temperature dependence of the resonant frequency of the fork in this tem- 
perature range comes from the viscous mass enhancement in Eq. 1)20)1 . Our 
experimental data for normal 3 He together with a fit to Eq. 1)20(1 are shown 
in Fig. 1101 The quantity plotted on the vertical scale is the term in the 
parentheses on the right hand side of Eq. (|20|) . The fit gives (3 = 0.90 and 
B = 0.73 for the mass enhancement parameters of the fork used at zero pres- 
sure and (3 = 1.25 and B = 1.05 for the fork used at 29 bar. The (3 factor 
can also be determined from measurements at the very lowest temperatures 
in 3 He-B, well in the ballistic regime, when the shift of the resonant fre- 
quency is caused entirely by inertial effects. This way we obtain (3 = 0.88 
for the fork used at zero pressure and (3 = 1.20 for the fork used at 29 bar. 
These second values are close to the ones shown as the zero intercepts on 
the vertical scale of Fig. 1101 

We conclude that the mass enhancement factors (3 and B turn out 
to have rather different values for our two forks, in spite of the fact that 
these two forks have closely similar dimensions and room temperature os- 
cillator properties. In both cases (3 is larger than the theoretical value 
(3 = (7r/4)W/T = 0.69 for a single beam of the size of one leg, as might 
be expected owing to the two legs of the fork in close vicinity of each other, 
Sec. 13.1.1 The sizeable differences in the fitted parameter values do not sup- 
port our hopes that quartz tuning forks could be used as reproducible sec- 
ondary thermometers, with a common calibration for all forks of the same 
type. 

Another important question is whether for a given fork the calibration 
obtained for one pressure can be used to interpret results at another pres- 
sure without re-calibration. To check this we repeated measurements at zero 
pressure using the fork for which a calibration at P = 29 bar had been ob- 
tained. In the regime 5 <C (T>,T,W), where the simple model from Sec. 3.1 
is applicable and which at zero pressure corresponds to T > 3T C , the mea- 
sured resonance frequency and width are within 10% from the prediction of 
the model. When temperatures approach T c the deviation increases and at 
T c the resonance width is twice larger than expected. In this temperature 
range evidently the interaction between the two legs of the fork becomes 
important and probably the surrounding of the fork (which here includes 
another fork less than 1 mm away) also influences the result. Thus we can 
conclude that the simple hydrodynamic model presented here describes rea- 
sonably well the behavior of the fork in normal 3 He but is not sufficient 
for exact scaling of a temperature calibration from one pressure to another, 
especially in the regime of large viscous penetration depth close to T c . 
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Fig. 11. Resonance frequencies of two quartz tuning forks in liquid 4 He at 
saturated vapor pressure, measured in two different setups. The symbols 
represent experimental data. The lines have been fit to the data at T > T\, 
using Eq. (|2U|) with (3 as fitting parameter and B = 1. The fit gives (3 = 1.39 
for setup 1 and j3 = 1.27 for setup 2. The data on the physical properties of 
liquid 4 He are from Ref. |26| . 

7. TUNING FORK IN 4 He 

Our measurements in liquid 4 He have also been performed in two differ- 
ent setups. In both cases the forks^ are partially inside their original cans, 
only a hole is ground in the can to provide a connection between the fork and 
the liquid in the 4 He sample container. The temperature is determined from 
the saturated vapor pressure. In the first setup the LHe temperature vac- 
uum parameters of the fork are /o V ac = 32708 Hz and A/ vac = 0.4 Hz. In the 
second setup a fork is used with /o V ac = 32709.97 Hz and A/ vac = 0.06 Hz. 

The resonance frequencies of the two forks are shown in Fig. as a 
function of temperature. In liquid 4 He above the superfluid transition the 
density changes faster than the viscosity. Thus the resonant frequency is 
more useful for thermometry. Indeed the measured resonant frequency is 
reminiscent of the inverse of the well-known liquid density, with a maximum 
in the density just above the superfluid transition T\. However, a fit of 
the measured frequency in the normal phase to Eq. (J20[) shows systematic 
differences which so far remain unexplained. In normal 4 He the viscous pen- 
etration depth 5 is of sub-micron size. Thus the influence of the small viscous 
term in the added mass cannot be reliably distinguished in the presence of 
the rapid variation of the liquid density. Therefore, in Fig. ^2 the value of 
B is fixed to 1 and the only fitting parameter is (3. 
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Fig. 12. Scaled resonance frequencies of two quartz tuning forks in liquid 
4 He at saturated vapor pressure. Plotted in this way both sets of data from 
Fig. ^2 coincide. The quantity on the vertical axis is [(/ovac//o) 2 — 1]/ 
which can be considered an "effective" value of the density p/p q . 



Below T\ we plot the extrapolation of Eq. (|2U|) in Fig. ^2 assuming a 
simple two-fluid-model interpretation. The interaction of the normal and su- 
perfluid fractions via the possible existence of quantized vortices is neglected. 
The inertial contribution to the effective mass is attributed to the whole fluid 
while the viscous contribution is assumed to originate only from the normal 
fluid component. Thus (/ vac//o) 2 = 1 + (3p/p q + BS/ (Vp q ) i]p n /irf , 
where p n is the density of the normal component. In view of the systematic 
difference between the fit and the measurements in the normal phase, the 
extrapolation to T < T\ is very reasonable. 

In Fig. ^1 we plot the same data once more with [(/ovac//o) 2 — on 
the vertical scale versus temperature. Now the two sets of data are reduced 
on the same temperature dependence, which is mainly that of the total 
liquid density p. The nice agreement in this plot lends support to the model 
expressed by Eq. (|20|) and to the extracted values of j3 for the two forks, 
derived by fitting their data separately to Eq. 1)20(1 at temperatures T > T\. 
Similar to the 3 He results in the previous section, we obtain rather different 
values for the (3 factors of the two forks (f3 = 1.39 and 1.27, Fig.lTT|). These f3 
values are larger than those from the 3 He measurements, presumably owing 
to the can around the fork. Interestingly in Fig. I12( the resonant frequency, 
when scaled with the relevant (5 value, looks promising as a calibration for 
thermometry in 4 He at T > 1.5 K. 

In Fig.^Jthe resonance widths of the two tuning forks are shown. Above 
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Fig. 13. Resonance widths of two quartz tuning forks in liquid 4 He at satu- 
rated vapor pressure. The symbols represent the experimental data and the 
line is a fit of the data from setup 1 at T > T\ to Eq. (|21|) using C as fitting 
parameter. The fit gives C = 0.90. The fit to the setup 2 data is almost 
identical to the one shown and gives C = 0.88. The data on the physical 
properties of liquid 4 He are from Ref. |26j . 



T\ the viscosity of liquid 4 He is not a strong function of temperature and a 
fit of the data to Eq. (|2"T|) . using C as a fitting parameter, works reasonably 
well. Thus it is interesting to compare the extrapolation of the fit to the 
experimental data in the temperature regime T <T\. In the extrapolation 
we use the same model of non-interacting normal and superfluid components. 
The resonance width is only associated with the normal component and thus 
in Eq. (|21|) we replace p with p n . Remarkably, the two forks follow nicely 
the extrapolated dependence in the superfluid regime, in one case down to 
1.8 K and in the other to 1.4 K. 

An intriguing question is whether any contribution in these results can 
be attributed to vortex generation and mutual frictionP^ From earlier mea- 
surements with vibrating wires, grids, and spheres it is known that extra 
damping occurs in superfluid 4 He even at low drive from the interaction of 
the normal and superfluid components in the presence of quantized vortices, 
from mutual friction losses. In 4 He vortices are easily pinned on surfaces as 
the vortex core is of atomic size and any surface becomes sufficiently rough 
for pinning. Thus the surface of the vibrating sensor might be loaded with 
pinned remnant vortices. These vortices might originate from thermal coun- 
terflow produced in a rapid cool down through T c , from some other source 
of residual flow in the system, or if the oscillating object has been driven 
previously in the superfluid state at velocities above some critical value. 
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In the present measurements shown in Fig. 1131 the amplitude of the 
fork current was always kept low (typically below 120 nA). This current 
corresponds to velocities less than 1.5 cm/s, which is below the typical critical 
velocity in 4 He of order 4 - 5 cm/s. Additional tests were performed when 
the forks were driven hard enough so that non-linear response could be seen 
and vortices should have been created. After the drive was reduced back 
to low level the width also returned to the original value shown in Fig. 
without hysteresis. This stability in the response is quite unlike the usually 
observed differences between "virgin" and "trained" responses of a vibrating 
grid, for instance.^ If this stability of the forks can be reliably reproduced 
and especially if it persists down to lower temperatures, then the quartz 
tuning fork might become a useful thermometer for superfluid 4 He which is 
not troubled by "vortex layer" problems. 

8. CONCLUSIONS 

The quartz tuning fork is a robust and easy-to-use sensor in cryogenic 
environments. In view of the results in Fig. [7J it appears to be a useful 
secondary thermometer for superfluid 3 He research. It requires less work 
and know how to implement and to operate than any other of the currently 
available methods in the superfluid 3 He temperature regime. Its response 
is well described at low excitation in the linear drive regime in terms of 
our hydrodynamic model which includes fitting parameters with a physical 
origin. More statistics on different forks are needed to decide whether simple 
means can be worked out to fix the parameter values and the calibration of 
the device as a thermometer. For such tests the forks should be adequately 
thermally cycled and preselected based on their LHe temperature resonance 
widths A/ vac in vacuum. The important physics, which should be explored 
with the quartz tuning fork, lies at the lowest temperatures in the ballistic 
regime, both in 4 He-II and 3 He-B, where the interaction of the fork with 
quantized vortices should be investigated. 

EPILOGUE: Dedication to Frank Pobell 

Oscillating devices immersed in a bath of liquid He were an important 
element in Frank Pobell's research. He is remembered for his passionate 
mission to reach ever lower temperature records. A controversial element in 
this quest was thermometry. The vibrating wire resonator was and still is the 
best thermometer for the T — > limit in 3 He-B. To explore the limits of this 
device, he studied in many papers the nonlinear response and dissipation in 
different wire materials down to below 1 mK in vacuum. He discovered that 
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even metals display slow heat release from defect structures, which relax 
similar to the tunneling model of two-level systems in glassy amorphous 
materials.^ These measurements should now be repeated for the quartz 
tuning fork. 

The senior members among the authors of this report remember Frank 
from this time as an extremely focused and industrious researcher. He was 
a visitor in the Low Temperature Laboratory for two months during the 
spring term of 1991. One of us (MK) was sharing the office room with 
him. He was working from early morning until late evening without break 
on his administrative chores, writing his research reports, and examining 
the manuscripts submitted for publication in JLTP. We were so impressed 
by this diligence and the steady flow of new results. An excellent example 
is his booJpQl "Matter and Methods at Low Temperatures", which he had 
just completed and which appeared later in the same year. This book is 
still gratefully used as the best text book for our courses in low temperature 
physics. East of Germany Frank is remembered (LS and PS) for his help and 
support to the low temperature community behind the iron curtain before 
it was lifted in 1989. This special relationship has continued over the years, 
an example was Frank's co-chairmanship of the LT21 Conference in 1996 in 
Prague. 
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